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ABSTRACT 

The design and operation of a special purpose digital signal 
processiiig system is explained. Tlie systesn consists of a signal 
analyzer, an interface between the analyzer ard the M.I.T. Compat- 
ible 'Time Sharing System, and appropriate computer progi'aras. With 
iniuxts from an instnur^nted shaft in the M.I.T. variable-pressure 
propel3,er tunnel, the system permits rapid spectral analysis of 
propeller generated forces and moments in six degrees of freedom. 

Thesis Supervisor: Justin E, llenrir 

Title: Professor of Naval Architectxure 
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1. INTRODUCTION 



1.1 Background. on July l, 1968 , the Department of Naval Archi- 
tectTore and Marine Engineering at M.I.T. sii.braitted a proposal to the 
National Science Foundation to build a system that would permit un- 
steady hydrodynannc measurements in the variable-pressure vrater tunnel 
operated by the department (l) . research into means of reducing vib- 
rator’y excitation to the hull of a ship depends upon the availability 
of methods to identd^y and measure periodic signals of Icnown frequency^ 
periodic signals of unlcnoim frequency, and random signals. Of parti- 
cular interest are the six components of periodic forces and raoments 
on the propeller shaft due to non-xmiform flow about the propeller and 
the fact that the propeller has a finite ntunber of blades. 

The design and construction of a specially instrumented shaft 
vjhich xrill enable simultaneous measurement of the instantaneous values 
of the three forces and three moments due to the propeller interactiori 
is covered in the Naval Engineer's thesis by Charles 0. Horton (2>. 
A matrix of miniature strain gauges and integrated circuit operational 
amplifiers built inside the shaft provides six signals in the 100 mv. 
to 1 volt range at an impedance of a fevj ohms. All the signals of 
interest vri.ll be below 1 khz. Each of the six forces and moments to 
be stud3.ed is a linear combination of these six signals, and the mat- 
rix of coefficients for this compxxtation is to be determined by dyna- 
mic calibration of the shaft. A separate circuit in the shaft housing 
delivers a sharp synclironizang pulse once per revolution. Fovjer to 



3 ' 









I 



the amplifiers and output signals from the ajralifiers bxq transferred 
between the shaft and the outside via low noise slip rings. 

The design of a simple digital data~processing system to rapidly 
analyze the signals received from the instrumented shaft is the s^^b•' 
ject of this I'eport, 

1.2 CHuICiS OF SISTM. The principal requirement for the data- 
processing system \rere wide flexibility in types of analysis available 
to the experiir^enter, short tui’n-aro\md timoj and cost below $1^,000. 

Originally it >/as felt that all signal-processing equipment 
should be self-contained and located in the propeller tunnel experi- 
rrental area. Complete signal-processing systems which could correlate, 
signal average, and spectrum analyze turned cat to be too expensive. 
Considered next was a method of rGCordD.ng data from the experinent, 
converting it to dig3.tal form, and then pxmcloing it on paper tape or 
cards for submittal to the H.I.T. Information Processing Center for 
"batch" processing, -The long turn-around tirae rrn.de this prospect dis- 
coirraging. The final choice was to use a commercial digital signal 
analyzer in the experimental area, and to malce an interface between 
the analyzer and the IHI 109h coir^mter used in the M. I. T. .Compatible 
Time Sharing System (CT3S) , The analyzer would collect data frorn the 
experiment, make the analog-to-digital conversion, and either store 
the data from a single sweep for analysis by an auto correlating or 
crosscorrelating prograra in the or malce hundreds of sweeps in 
averaging mode, leaving in its laemory a digitized signal of enhanced 



6 







I 



.'S-in 



( 



signal“to**noise ratio for use in a digital Foiu’ier analysis. 

'itie signal analyzer chosen irzs the northern Scientific Model 
with optional specifications as follows: 

Memory: 102li. addresses; 10^ coimts per address 

Number of jjaputs; 8 (with NS~30p eiglit input expander) 

Readout: Te3-etype Model 33 ^vSkTC (includes tape 

pimch and reader) 

The met the requirements of this project in the follovring 

specific ways : 



1, Becavise it could read its output into a Teletype;) which in tiun 
can be coimected to CxSS i-7ith a siiqile electrical DUiterface to a Data- 
Rhone, the cox^^utcr interfacing problem would not be coifrolj.cated. 

2o The system, in the irorst case, should resolve the thii'd harmonic 
of a 9~bladed propeller operating at 2000 r.p.ra, with all 8 inputs 
similtaneously. With time muH-tiplexing of the eight inputs, this 
requires a sampling frequency of at least llu^t khz,, or . a-.: sampliing 
period of no grea/ber than 70 microseconds. The NS^^^O can sai^le at 
intervals of 50 microseconds. 

3. The total price of the analyzer and peripheral gear vias aboiit 

$ 12 , 000 . 

The electrical intei'face betvjsen the Teletype machine and the 
Da'ba-Phone vjas designed and built by the author. It is electronically 
very simple and requires no attention or adjustment. 
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Block Diagram of the Data Processing System 
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The system is coriipleted by a conversational prograju written in 
Fortran TJ which receives typed input froBi the experimenter and the 
signal analyzer and returns a Fourier analysis of each of the six 
forces and moments. For an analysis involving all 3,0 2l| memory loca^ 
tions in the analyzer, the turn-around time is about 20 minutes, most 
of which is aaialyzer read-out time. The program listing and docxirren- 
tation is given in Appendix A, 

J?'ig, 1 is a blocJc diagram of the completed signal analyzing 
system. 
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2. MODIFIGATiONS TO EQUMlEiNT 



2.1 INPUr MODIFICATIONS. The Na~550 analyzer has an input iirpadance 
of at least lOOK and sTri-tch-selected AC or DC coupling. The torque 
and thrust outpvits from the shaft i-ri.ll have d.c, levels considerably 
higher than tiri.ce the a.c. signal components ^ and consequently the 
baseline adjustments on the analyzer vrill be insufficient to remove 
them. Unfortunately, the time constant of the a.c. coupling is .01 
second, v;hlch is too small to pass low-frequency signals from the ex- 
periment without significant attenuation and phase shift. ^ This con- 
dition can be corrected by replacing the .1 pf, coupling capacitors in 
the eight input circuits vrith 10 jif . capacitors of equal voltage ra- 
tings . 

2.2 SI'JSBp modifications. VJith 8 inputs and a memory cycD.e time of 
50 microseconds, it takes ,U msec, to add one data point to each of 
the eight memory subdivisions. For simplicity, assume that each chan- 
nel uses 12^ memory locations. Vie desire to store data representing 
one complete propeller shaft revolution in each subdivision. The 
maximum period of one I'evolution .is then 12^ X .i^'msec,, which corres- 
ponds to a minimum shaft speed of 1200 r.p.m. Increasing the speed 
does nothing more than reduce the number of data points corresponding 
to one period. 

Thus, for speeds of 1200 r.p.m. and above, the time base of the 
NS-5^0 must be set for .1| ms, /point. Similarly, for speeds betxjeen 
600 and 1200 r.p.m., a sett.ing of ,8 .liis./point is required. However, 
as siipplied by the manufacturer, the NS -550 skips over this range. 
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going from ,2 ms. /point to 2 ms. /point. 

The required modifications to the time base circuitry are not 

difficulty and circuit diagi’ariis shoifing these changes can be obtained 

2 

from Northern Scientific. 

2.3 LINE FEED HODIFICATIONo . V/henever CT33 is called iip from a Tele- 
type terminal, the corputer begins by signaling the automatic anst-rer- 
back device on the Teletype. The codes which are returned identify 
the particu-leir Teletype machine and tell the con^Duter if it must sup- 
ply separate line feed signals whenever a carriage return is either 
sent or received. The two codes are identical except that in one a 
letter "G" appears in the two places which contain tlie letter "D" in 
the other code. If the ans^xer-back device is set up for the "D" code, 
the computer supplies a separate lino feed. If it is set up for the 
■"C" code, the coinputer does not supply the separate line feed, but as- 
s-umes that whenever a carriage I’eturn is either sent or received, the 
Teletype machine will line feed on its ovxn. 

The "D" code works satisfactorily for all input and output ex- 
cept for signal analyzer readout, because the signal analyzer produces 
a line feed of its oim, and does not wait for the receipt of the com- 
puter-generated Line feed before continuing. This causes confi\sion 
and loss of data. 

Using the "C“ code solves the above problem and introduces ano- 
ther: the operator must type his o^m carriage retm’ns when typing 
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input to the computer, and must irjanually advance the paper for each 
successive line of data returned by the computer. 

The solution to this problem is to set up the fvmctlon box in the 
Teletype machine for "Suppress Line Feed" and "Neu Line, Lino Feed," 
in accordance tTxth the technical manual supplied mth the machine (3j. 
The "Suppress Line Feed" feature ignores all line feed signals re- 
ceived, and the "Keir Line, Line i-eed" function automatically produ- 
ces a line feed vrhenever a carriage return signal is either transmt- 
ted or received, Uf course, the “C" answer-back code is used. 
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3. electrical interface BEE-ffiEN TELETIPE AlID BaTa-PHONE 

3.1 TELETjTi"E SIGMLS, xhe 33ASRTG Teletype that Northern Scientific 
supplies as an optional readout device for the NS“550 coiaes wired by 
them for ttiree modes of operations “Tape to Typej," "Analyzer Read 
In," and "Analyzer Read Outi" In the first mode, the tape reader sup- 
plies data to the printer. In the second mode, the tape render sup- 
plies data only to the memory of the analyzer. The D.ast mode has the 
printer (and tape punch, if desii’ed) print (punch j o\i.t the data in the 
analyzer’s nsmory, . In all these m.odes, the data is transferred in 
pai'allel to and from the analyzer. The parallel- to-serial conversion 
of the code for each symbol (required for transmission via the Data- 
Phone) is acconplished electromechanicaU.y by the printer mechanism. 

The Teletype itsei-f has px’ovision for both "local" and "line" 
operation. In local m.ode, an internal power supply provides the 6o 
ma, loop current required to operate the selector magnet system. When- 
ever the machine is being used for purposes not requiring a connection 
to the coTipixter, it shou3.d be operated in local mode. In line mode, 
an external povrer supply is required, and it is connected to terminals 
6 and 7 of terminal strip X, (For details of the Teletype wiring, see 
reference (I|.) . If the current from the external supply is modu-lated 
on and off by the ASCII code f oi' a letter , the mach3.ne will print the 
letter (and punch it as ire 3.1, if desired). Alternatively, i.f a key on 
the keyboai'd is ■deppressed, the Teletj’pe itself modulates the line 
current with the ASCII code. Thj.s also happens irhen a symbol is 
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Fig. 2 



Sch^natic d.isLg^L'OJn of the electrical interface 
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printed during operation in “Tape to Type" or in "Analyzer Read Out" 
mode . 



3.2 DATA-rHOI'IE SIGNALS , The Data-Phone transmits and receives all 
data serially, in the form of pxilses of (norainal3.y) pins and minus 3.5 
volts. Full details are given in tlie manufacturer's interface speci- 
fication (5^ . 

3.3 THE INTERFACE . The interface translate.? the modulated ]jLne cur- 
rent from the Teletype to positive and negative volta-ges for the Data- 
Phone, ajid vice versa. Fig, 2 is a schematic diagram of the unit. 
The interface consists of a 6o m.a, loop supp3.y, a ±10 volt supply, and 
tu-o transistorized stTitches, The incoriiing signal appears at terminal 
3 of the "customer's equipment" connector, Mienever the "Hai'k" (nega- 
tive) signal, is received, the loop current is svritched on by Ql, - The 
"Space" (positive) signal causes Ql to svritch off. The germanium 
diode connected to the base of Q2 serves tuo pxu^poses. First, its 
back resistance is sufficient to asswe that Q2 is cut o.ff whenever 
the Data-Phone is disconnected. Second, it prevents damage to the 
base -to-emit ter junction in the event the Data-Phone produces a nega- 
tive voltage greater than -10 at terminal 3. (It may go as high as 
-25> according to the interface specification,) 

Q3 and Q!; supply the outgoing signal. The vo3.tage divided con- 
nected to the bases of Q3 and Ql} is chosen so that when the Teletype 
is marking (terminals 6 and 7 of strip X shorbed) a negative voltage 
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at loTr Djupedance is presented at pin 2j and v;hen the Teletype is spa~ 
cing, a corresponding positive voltage atplovx jjnpodance appears ..at ^pin 

2 . 

The continuous positive voltage connected to pin 20 enables the 
Data-Phone to enter and remain in the "Data" riiode. 
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U. PROGRAI-MING REQUIREMENTS 



li.l OUTPUT FORMAT FOR THE NS-<^0. A typiccal line of output from the 
NS-5^0 consists of a line feed (LF), a cai’riage I'eturn (CR), and eight 
6-digit numbers, each separated by a single space from the preceeding. 
If each memory location contains only decades (which is the case for 
the NS-^50 used in these experiments), the first digit of each number 
is always zero. The analyzer never prints a plus or rdinus sign, V/ith 
the modification described in section 2.3 of this report, the Teletype 
ignores the line feed, but automatically produces one when it receives 
the carriage return. 

The initial OR and LF is interpreted by the computer as a blank 
card image, and in order to avoid having the first 8 locations of the 
input array set equal to zero, the program mu.st be vrritten to read the 
blank iiiiage before reading the first 8 data points. 

The following program segment vjill successfully read JMAX points 
from the signal analyzer: 

READ 100 

100 FORMAT (IX) 

DO 10 N=l,J>kX, 8 

10 READ 101, (lNDATA(NZ),NZsN,NA) 

101 FORi'IAT (8(I6,LX)) 



17 



14.2 AjJLOCATION of storage in the NS-550. The 102)4 raeiTiory locations 
5 ji the analyzer are divided equally among the input channels being 
used. The first and last locations in each block of raemory are alirays 
set to zero^ except for the first location in the first block, which 
always contains the number of sweeps made in the experiment. uata is 
always read out sequentially by memory address number. Data from Chan- 
nel A appears in the first block of. memory, data from Channel B ap- 
pears in the second block, and so on, A computed GO TO statement vrith 
its argum.ent equal to the number of input channels can be used to as- 
sign the data received from the analyzer readout to arrays indexed by 
channel number (See Appendix A) . 

14.3 INPUTS aUrPLIED BY THE EXPERH'jENTER . Becaxise tDrds -Sharing allows 
conversational interaction between the program and the experimenter, 
the prograjii can ask for the values of various constants and receive 
the typed responses directly. The follox-iing sets of constants are re- 
quired: 

1, The number of input channels used in the experiment 

2, The settings of the input attenuators on the analyzer 

3 , The matrix of coefficients relating the actual forces and moments 

3 

to the signals appearing at the input of the analyzer 
I4. If the experiment involved signal averaging, and the readout did 
not begin vjith the first address of rcem*ory: the nximber of sweeps 

made by the analyzer 

If the program does a Fourier analysis; the number of data points 
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in a channel which represent one coriip3,ete shaft I'evolutionj and 
the number of harmonics to be computed 

6, If the program compubes a correlation fimotion; the range of the 
delay time (the arguraent of the correlation fraction) 

7. If the program crosscorrelatcs: the values of the second fraction 

or :• an algorithm for computing these values 

The Fourier analyzing program in Appendix A has nuEierous examp3,es 
of conversational jnput-output routines. 
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OPERA.TIKG INSTRUCTIONS FOR A SIX»iJll®ISIONAL AMLYoIS 



It is assumed in these instractlons that the reader is famliar 
with Ilr, PIca'ton's thesis (2) which describes the overall design of the 
complete experimental system, the operation of the NS- 550 , which is 
covered in the instruction manual (6), and the CTSS procedures, which 
are covered in the CTSS User's Manual (7). 

PROCEDURE 

5.1 Push the STOP button on the NS-550, 

5.2 Set: DECADE SELECT to 10^ 

READ MODE to CRT 
aDD~SUB to ADD 
MEI'IORI CROUP to 1/1 
TII'IE BASE to 2 ms. /point 
NO. OF SIEEPS TQo^^ 

DEMY TD-IE to 1 ms. 

FRE DEMY,' POST DEMY to OFF 
ADDRESS SELECT to OFF 
VERTICAL OFFSET to OFF ' 

COWERSION GAIN tv nVE 

SVffiEP MODE to AVE 

TRIGGER to RECURRENT 

VIEW IWPUT-Va.EW MEMORY to VIEW MEMORY 

INPUT SELECT to nXT 
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5.3 Push START READOUT, a pattern of some type should appear on the 
CRT, If nothing appears , carefully adjust the INTENSITY, HORI- 
ZONTAL, and VERTICAL controls tuitil somethrlng does, Al.ways keep 
the INTENSITY control set to the minimum level which produces a 
visible trace or spot on the CRT, Tnis is especially important 
when only a single spot appears on the CRT. 

Push the ERASjfi button. The CRT pattern should imimediately degen- 
erate to a straight line. Adjust the VERTICAL POSITION control 
so this 3.ine appears at the very bottom edge of the CRT face. 
Adjust the HORIZONTAL POSIT.i.ON and HORIZONTAL EXPAND controls so 
that this line extends sDdghtly less than the vridth of the CRT 
face , 

5.5 Push STOP. Tnen piish START MEASURE. Let the iinj.t run for about 
3.0 seconds, then push STOP, followed by START READOUT, Tliis cas- 
ses a small amount of noise to be read into the memory. 

5.6 Set the DECADE SELECT to 10^ and observe that 10 horizontal lines 
appear on the CRT. Using the VERTICAL EXPAND and VERTICAL POSI- 
TION controls , male these 10 lines coincide with 3.0 o.f the 11 
scribed reticle lines. This coiiraletes calibration of the CRT. 

5.7 Push STOP. Set DECADE SELECT to 10^. Set the INPUr ATTENUATORS 
for all U inputs of the NS-550 'and all k inputs of the NS-305 
to 1 volt. Set all 8 of the AC-DC siri.tches to AC. 

5.8 Connect the six stra3ri-gauge outputs from the shaft to channels 
A tlxrough F. Connect the synchronizing pu3.se cable from the 
shaft to the positive input of Channel H and also to the EXT 
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TRIGGER connector 



^,9 Start up the props llei' and set the shaft r.peiUc to the value that 
ull3. be used in the exj^sriment . Set the TEE BaSE to ms,/ 
point if the r.p.ra. is above 1200, Set the TEE BA.SE to ,8 ns./ 
point if the speed is betvraen 600 and 1200 r,p.m. Set the TEE 
BASE of the NS-305 to the same value that is set on the NS=-550, 
$.10 Set VIE’.-/ EJPUT to A. Set TRIGGER to -{'MT, and push START MEASURE. 
Adjust the TRIGGER SENSITE7ITI fox’ stable triggering of the hori- 
zontal sweep i as indicated by a steady pattern on the CRT that 
does not raove left or right. Center the pattern vertically by 
xxsing the Chamiel A BASnLEE ADJUST, If the pattern goes off the 
top or bottom of the CRT, set the Channel A E^rUT ATTEl'JUATOR to 
10 volts. If the pattern is very small, set the Channel A EIPUT 

ATTEMATuR to 0,1 volt. 

^.11 Sxritch VIEW E^rUT to B. Using the Channel B BASELDE ADJUST and 
INPUT ATTENUATOR, position and vertically scale the pattern as 
before, nepeat for Chamxels C and D. 

^,12 Set the DEC,J)E SELECT to LOG, and VIEW INPUT-VIEW 1-M-IORI to VIEiJ 
MKI-IObY (xlxis periTiits viewing of the foiu’ inputs of the NS-30p). 
Set VIEW INPUT on the Mo-305 to E. using the Channel E INPUT AT- 
TENUATOR and baseline ADJUST, position and scale the pattern as 
before. Repeat for Chaime3. F, 

^.13 Bet the VIEW ELUT to G, Center the line using the channel G 
BASELINE ADJUST. 

Set VIE.'/ INPUT to II, The pattern should be a straight horizonl^l 
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line except for a pvJ.se near the right-hand side. V/ith the EA.SJ:;- 
LDfB adjvist for Channel H; position the straight line about IJj^i 
up fj-’orn the bottom. With the IWPUi ATiEhTJATORj make the pulse 
height at least two divisions on the reticle, 

^.1^ Push STOP. Set DECADE SELECT to 10^, Push START PlMDOUT. Pv.sh 
ERASE, rush STOP. This completes the adjustments, 

^.l6 To begin the measurement, push START 1-IEASURE, Initially the CRT 
will shovr only a horizontal line at zero. As the experiiaent con- 
tinues, vertical deflections will apperir and grow. Set DECADE 
< 

SELECT to 10 . Continue the e:xporiment until the position of any 
portion of the trace approaches the top of the CRT, Then push 
STOP, (Alternative; Set the nO. OF S1JEEP3 control to 800 before 
pushing START MEASURE, The analyzer will msQce 800 sweeps and 
stop automatically.) This con^letes the raeasurement phase of the 
experiment , 

5.17 Turvi on the electrical interface between the Teletj’pe and the 
Data-Phone, Turn on the Teletype in "Line“ mode. Set the func- 
tion m-jitch to "AMALYZER READ OUT." Dial up the computer and 
LOGIN. Then type LOADGO PROP. The program Trill begin by asking 
for a nuraber of constants. The .follovring steps give the details. 

5.18 "TYPE DIUl'IBER OF INrUT CHANNELS." Format is II. Type in the num- 
ber six. (The sync input on Channel H is not. counted as an in- 
put.) 

5.19 "type INPUT ATTENHiTOn SETTINGS." Fonnat is 6F7.2 Type in the 
settings of the six INPUT aTTEIRJaTORS , A through F. 
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Foriiiat is II. Type 6, 



^,20 "TY.E NUl-IEESR OF OUTPUT CmillELS DSSIRk;D." 

5.21 "TTPS UUl-IBiiR OF HAK-'ONICS DESIRED." Format is 12. The basic 

period of analysis is one shaft revolution, EXAl'ir'LE: To obtain 

analysis up throuigh the third hariaonic of the blade rate for a 5“ 
bladed propeller, 15 har-monics are required. The maximun nuniber 
that can be typed in is 50. 

5.22 "EKTSR TrIE M-OUT BY N-IN CALIBRATIOII MATRIX." The response must 
be the 6X6 strain-gauge calibration matrix , where k, the 
row subscript, is the output channel number, and i, the column 



subscript, is the input channel niuabs 



The dimensions of each 



entry must be either ft ,-lbf ,/volt or Ibf./volt, and the conputer 
ijill read s?sc sixccessive rot-rs of six numbers each on a 6 f 7.2 for- 
mat. 



5.23 "TYfS ITOHESR of D:ATA points to EE READ." Forxaat is ll|. Type in 
the number 10 2U. 

5.2li "HERE IS HOW I HATO YOUR SETTOTGS-" The coirputer vri.ll read out 
all the values as it has them, and gives self-explanatory in- 
strxictions on what to do if there is an ei'ror anyvrhere. 



NOTE: In the next two instructions, the program mentions svri.tch- 
ing to eitlier mamxal or automatic line feed. If the modifica- 
tions described in section 2,3 of this report have been made to 
the Teletype, the line feed instructions should be deleted from 
the program, 

5.25 "I M PJLADY TO ACCEPT DATA FROM THE SIGNAL AVERAGER." If a paper 



tape .of the readout is dSvSirod, turn on the tape punch and hit 
the space bar on the keyboard about 10 times to raake a leader on 
the tape, but do not type a carriage return or a line feed. Push 
STOP on the 1IS“^^?0, Set the READ MODS to TIPS and push START 
REivDOUT. The readout Trill take about l5 minutes. 

^.26 "SI-/ITCH TO AUTOI'iVTIG LINE FEED AND TYPE GAPJlIAGn RETURN." See 
the note above. If this instruction and the read statement im- 
raediately follorring it in the program (Appendix A) have not been 
deleted, sirrply hit the carriage return to continue, 

^.27 “TYPE NUMBER OP D'lTA POINTS COMPRISING OnE CuMPLSrE PERIOD." In 
the event the program was not able to detect the sync pulse in 
Channel 8, it trill rna]ce the above request. If tliis happens, look 
at the last 16 ro>JS of the printed readout. Observe that the 
sixteenth row from the bottom begins vrith 000000. Number these 
rows 1 throitgh l6 starting vrith the one begjjnning with zeroes. 
Then find the discontinuity in the printed data values. If R is 
the rovr and C is the colvunn viiere the discontinuity appears, the 
number to type is 8(R»ljd-(G-l) , The format is lii. 

^,28 "REaD uUT NUi'iBER OF PnSSES." If the computer receives a zero for 
the first numbex’ read in from the analyzer, it vrill make the 
above request. Push STOP on the Mo-5^0, Set 1M40RY GROUP to 1/1 
and turn the ADDRESS SELECT svritch OFP. Set RE/i.D MODS to 'CRT, 
Push START READOUT and then push STuP. Set READ MODE to TYPE and 
push START REilDOUT. After the first nurabei* is printed and before 
the eighth is printed, push STOP, set READ MODE to CRT, and then 
hit the carriage return. 
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( 

$.2? Tlie compu.ter ixill now read out by channel number and haj'monic 
number tho amplitude in Ibf . or ft»~lbf and the phase dn deg=- 
rees of each harmonic. The zeroth harmonic is the average ^ and 
it has no significance in this type of experiment, since the in^ 
puts are all a.c, coupled. The channel nurabsrs can be related to 
a specific force oi’ moroent by a brief consideration of the compuc^ 
tation inplied by the 6 A 6 calibration matrix, 

5.30 "THIS COIIPLETES THE AMLYSIS. ITPE 1 TO START AU-IBI, 0 TO STOP." 
If a "1" is typed, the prograi?i nill return to step 5*25. 
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APPEl'IDIX A: "PROPj" A GOm-'Ul'BR PROGRiJi FOR FOUitlER AMLYSIS 

The program described in thj.s section embodies the principles 
outlined in section h s,nd is the progrejn used for the 6“d±m3nsional 
analysis of section 5. The listing begins on page 32. PROP can be 
divided into ^ pa.rts: 

1, TYPED INPOTS: PROP requests values from the experimenter and pro- 
vides an opportunity to start again if errors have occurx’ed. (Cards 
1-37) 

2, AHaLYZER READOUT: The program reads the analyzer output into the 
array IIRDATA^ and then assigns the values to appropriate blocks in tiie 
array UJPUT, according to the rule stated in section U.2. If 3' or 
more inpu.t charnels are being used, the program scans the contents of 
Channel 8, looking for a discontinuity of at least 20,000 counts. If 
it finds such a discontinuity, it sets the number of data poltits rep- 
resenting one complete shaft revolution equal to the location of the 
discontinuity. If thb discontinuity is not fo^^nd, or if the number of 
inputs is or less, PROP i*equ8sts the experimenter to type in the 
correct value, (Cards 38"87) 

3, MATRIX MULTIPLICATION: The program converts strailn-gauge signals 
into forces and moments. (Cardjs 88-S>8) 

It, FOURIER AMALYSIS: PROP continues by calculating the coefficients 
for the Fourier sine and cosine series for each force and moment, 
(Cards 99”H6) 
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OUTi^UT: Th3 program determines the number of signal analyzer pas- 

ses ( sweeps )j calculates the magnitude and phase of each harmonic for 
each output channel, prints the results, and provides the experimenter 
the opportxmity to read in neu da.ta for the nex"b analysis. (Cards 117- 
152 ) 

The folloT-Ting is a chart of the var*iables used an tne program. 



VARIABLE MAIIS 


RANGE OF values. 


DSSCRIrTIOh 


A(K,L) 


Ibf./lOO counts or 
ft.-lbf./lOO coxmts 


Fourier sine coefficient for 
the Lth harmonic in the Kth 
output channel 


Al-iP 


ft.-lbf. or Ibf. 


The magnitude of a harmordc 


ave'"' 


ft.-lbf. or Ibf. 


The average value of the 
signal 5n an output channel 


b(k,l) 


Ibf./lOO counts or 


Fovirier cosine coefficient 




ft.-lbf./lOO counts 


for the Lth harmonic in the 
Kth output channel 


c 


27T/Bi'!DFC 


Coefficient used in Fourier 
series calculations 


calib(k,i) 


Ibf ./volt or 
ft.-lbf ./volt 


The k,i entry in the strain- 
gauge calibration matrix 


FuXSCA(l) 


.025-’100, volts per 
100 counts 


Input attenuator setting of 
input channel i 


I 


1-8 


Index giving iiput channel 
nixmber 


n-iAX 


1-8 


Number of input channels be- 



ing used in the experiment 
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II®ATA(N) 



IwPUT(l,J) 



0-99999 



0-99999 



J, JP, oPO, etc* 1-1022 



m\x 

K 

L 

M 



M 

NDPC 



NOUTCH 

iilj I'lY 
NZ 

OUTPUT(K,J) 



PASSES 



PHASE' 



l-102lt 

1-8 

1-^0 

i-5o' 

l6^ 2k i etc, 
1-1022 



1-8 

0-1 

1-102U 

Ib./lOO counts or 
ft.-lbf./lOO counts 

-99999, 



-90, to -5-270, degrees 



i‘he Nth data poant read in 
froin the signal analyzer 

The Jth data point in input 
channel i 

Index giving data point nxua- 
ter 

Nuinber of data points read 

Output channel index number 

Index of harmonic being cal- 
culated 

Numter of harmonics the pro- 
gram is to calcixlate for 
each output channel 

Index vjsed in reading IHDATA 

Number of data points per 
channel xrhicb represent one 
co 2 iplete period of the input 
signal 

Number of output channels 

J3ranching variables 

Index u.sed in reading INDATA 

The Jth data point in output 
channel k 

The number of sweeps made by 
the signal analyzer in ati 
e^qjeriment 

The phase of a harmonic in 
Inference to a sine function 
of the same frequency 
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RIH 


0 ,- 99999 * 


Real value of ilJPUT(l,J) 


RJ 


1.-1022. 


Real value of J 


SC0S(J)‘'““'"' 


-1,0 to 1,0 


Stored value of a cosine 




-1.0 to 1.0 


Stored value of a sine 



'“ME is printed as the 0th harmonic in the output format, 

’Joe 

PHA.se is initially calcixlated in rad).an3 from the arctangent lib- 
rary routine j and then is converted to degrees , taking account of 
the quadrant, 

"''“'SSIN and SCOS are used to avoid repeated computatioiis of the sine 
and cosine valites needed for the i‘‘ourier analysis of a gi.ven hai'- 
monic in different outpu.t cliaiinels. 
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APPEIIDIX B: KOTES 



1. Representing the low-in^Tedance soxirce as an ideal voltage sour ce j 
the transfer fnnction for the capacitor coupled input is 






flTrRC-f 



I -i-AwRC-f 

A straightfon-xard calcivlation shox-rs that the phase of H(f) is 

^ = tan'' ("zif ^/' ) 



and the mgaitude is 

I 

(2-TTRC/)^ 

With the shiai’t r.pem, set to 9^0^ f is l6 hz, and rxith RC .01 sec,, 
Zr^RC-f 1. , ^5' 




Increasing the value of RC to 1 second permits reducirig the r.p.m, to 
on3.y 96, Under these conditions, the phase shift is only 5,7 degrees 
and the magnitxido is dorm by a factor of only ,005, 



2, Northern Scientific, Inc,, 2551 Nest BeltDdxe, P,0. Box 66, Md~ 
dle'ton, Wisconsin, 535^2, Phone (608} 836"65ll. 



3, In the compixter program of Appendix A, each enti’y iia the calibra- 
tion riiatrix is assumed to bo independent of freqxxency. If dynamic 
calibration of the sha^t shoi-xs tliat this assuxqxtion is not valild, the 
program xdll have to be rexiritten to do a Fourier axialysis of the in- 
coming data before doing the conversion of strain-gaxxge signals to 
force and moment outputs, xhen each harmonic xiill have to be convert- 
ed separately i.n the program, using va3.xies for the calibration matrix 

entries that are appropriate for the frequency of the haimioni.c. 
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